QCD-Inspired Spectra from Blue's Functions 



Maciej A. Nowak 1 , Gabor Papp 2 and Ismail Zahed 3 

1 GSI, Plankstr. 1, D-64291 Darmstadt & Institut fur Kernphysik, TH Darmstadt, Germany & 
Department of Physics, Jagellonian University, 30-059 Krakow, Poland; 

2 GSI, Plankstr. 1, D-64-291 Darmstadt, Germany & 
Institute for Theoretical Physics, Eotvos University, Budapest, Hungary 
3 Department of Physics, SUNY, Stony Brook, New York 11794, USA. 
(February 1, 2008) 



We use the law of addition in random matrix theory to 
analyze the spectral distributions of a variety of chiral ran- 
dom matrix models as inspired from QCD whether through 
symmetries or models. In terms of the Blue's functions re- 
cently discussed by Zee, we show that most of the spectral 
distributions in the macroscopic limit and the quenched ap- 
proximation, follow algebraically from the discontinuity of a 
pertinent solution to a cubic (Cardano) or a quartic (Ferrari) 
equation. We use the end-point equation of the energy spec- 
tra in chiral random matrix models to argue for novel phase 
structures, in which the Dirac density of states plays the role 
of an order parameter. 



1. Lattice simulations have established that QCD 
undergoes a phase transition at finite temperature [Q. 
The precise nature and character of the transition are still 
debated. The simulations, however, show that the QCD 
phase diagram with three flavors is rather involved, with 
a strong dependence on the value of the current masses. 

Important insights into the role played by chiral sym- 
metry may be found in the quark density of states at 
zero virtuality. Banks and Casher Q have shown that at 
this point the density of quark states is directly propor- 
tional to the quark condensate in the chiral limit. The 
spontaneous breaking of chiral symmetry is followed by a 
structural change in the quark eigenvalues near zero vir- 
tuality. It was also observed that some generic spectral 
oscillations take place in this region j^J9], reminiscent of 
Airy oscillations [[|. 

An important question in QCD concerns the struc- 
ture and character of the Dirac spectrum during a phase 
transition, as induced by changing the number of fla- 
vors, the character of the color representation, the quark 
masses, the temperature or the baryon density, to cite a 
few. Since chiral symmetry breaking is encoded in the 
eigenvalue distribution at zero virtuality, it is clear that 
a chiral phase transition would affect quantitatively this 
distribution. Since massive QCD does not have, strictly 
speaking, a good order parameter, it is natural to study 
the behavior of the spectral density near zero virtuality 
as a potential alternative. The possible change in the 
quark wavefunctions at zero virtuality is reminiscent of 
a phase transition from a delocalized and hence coherent 
phase, to a localized and hence incoherent phase. The co- 



herent phase will be identified with the Goldstone phase 
in our case. 

An efficient way to get the Dirac spectrum has been 
established through the use of chiral random matrix mod- 
els [|||^] . Chiral random matrix models are 0-dimensional 
field theories, inspired by the flavor and spin structure 
of QCD motivated Dirac operators. They share much 
in common with "two- level" Nambu-Jona-Lasinio models 
0. By restricting the discussion to only constant modes, 
the role of symmetry and the lore of randomness become 
dominant and transparent. Results from random matrix 
models both in the microscopic |8j and macroscopic Q] 
limit have compared favorably with some lattice simula- 
tions of the QCD spectra pjJlOfl . 

The purpose of this letter is to use newly developed 
methods by Zee jllj regarding the addition laws of ran- 
dom matrix theory, to discuss the generic character of the 
Dirac spectrum in the macroscopic limit and quenched 
approximation. For a large class of chiral random matrix 
models as inspired by QCD either in the continuum or 
on the lattice, and instanton models, we show that the 
resulting spectra follow simply from the general lore of 
random plus deterministic. In section 2, we introduce 
the concept of Blue's functions and review Zee's argu- 
ment for the deterministic plus random (Gaussian) case. 
In section 3, we show how the Blue's functions can be 
used to condition the end-points of the spectral distri- 
butions, resulting into quadratic equations for positive 
end-points (Cardano) or cubic equations for positive end- 
points (Ferrari). In section 4 we apply these results to 
massive Wilson fermions in the one- flavor approximation, 
as well as massless QCD on the cylinder. The resulting 
spectra and singularities belong to the Cardano class. 
We show that the spectral distribution of QCD with 2+1 
flavors in the Wilson formulation, follows from a linear 
superposition of the solutions to Cardano class. Instan- 
ton models with off-diagonal order fall into the Ferrari 
class. Our conclusions and recommendations are sum- 
marized in section 5. 

2. The spectrum of a single random matrix was deter- 
mined long ago by Wigner and others ]l2]j . Generically, 
the averaged level distribution of an N x N random ma- 
trix H (Hamiltonian) with a probability distribution 

P(H) = ±eM-NTrV(H)} (1) 
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follows from the discontinuity of the resolvent 
v(X) = lim lmG(z = A + ie) 

7T £-^0 



with 



(2) 



(3) 



The average < ... > is short for averaging with the 
weight (|l|). For polynomial potentials the explicit form 
of (|^) is known |l3j. In the simplest case of a quadratic 
polynomial, the random ensemble is known as Gaussian 
with the Green's function 



G(z) 

The discontinuity of 
Gaussian ensembles, 

i/(A) = 



— {z — i\J 4 — z 2 ) 



(4) 



is just Wigner's semi-circle for 



X 2 



(5) 



The problem of adding random matrices can be re- 
duced by using Blue's functions, as recently discussed by 
Zee The Blue's function B(z) is just the functional 
inverse of the Green's function G{z). Operationally 



B(G(z)) = G(B(z)) = z 



(0) 



The Blue's function for a constant matrix c (random ma- 
trix of infinite weight) follows from the corresponding 
Green's function 
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G(z) = 

through the substitution z — > B(z) in (^|) 
B(z) =c+- 



(7) 



(8) 



Similarly, the Blue's function for the Gaussian ensemble 
with the Green's function (^) follows from (Q) in the form 



B(z) = z + - 



(9) 



If B 1 (z) and B 2 {z) are the Blue's functions of two ran- 
dom matrices, then the Blue's function B\ +2 (z) for the 
sum follows from the addition law [fll 



B 1+2 (z) =B 1 (z)+B 2 (z) 



(10) 



Hence, the problem of finding the spectral distribution of 
the sum of two random matrices boils down to the follow- 
ing algorithm : First, construct the appropriate Green's 
functions G\ and G 2 . Second, find their functional in- 
verses, i.e. Blue's functions B\ and B 2 . Third, construct 



Bi +2 using (|Io|). Fourth, functionally invert B\ +2 us- 
ing the definition of the Blue's function, leading to G\ +2 . 
The spectral density of the sum is the discontinuity of 
G1+2 along the real axis. 

To illustrate these points, consider with Zee |ll[| the 
case where H is the sum of a deterministic matrix Hd 
with eigenvalues e$, (i = 1, . . . , TV) and a random Gaus- 
sian matrix Hr. The Green's function for the determin- 
istic matrix Ho reads 



G D {z) 



1 



\ ^ 



N ^ 



(11) 



Through the replacement z —> B{z) the Blue's function 
B D is 



N ^ 



1 



(12) 



N ^ B D (z)-e l 

The addition rule gives 

B(z) = B D {z) + B R (z) - - = B D (z) + z (13) 

where the second equality comes after we have used the 
known form of the Blue's function for the random Gaus- 
sian ensemble (||). Substituting z — » G(z) in ( |l2| ) and 
( |l3| ) we get a pair of equations 
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N ^ B D (G(zj) - e, 



z = B D (G{z)) + G(z) 
from which we eliminate Bp (G)[] to get 



G(z) 



1 N 

-Y 



Z - G(z) - 6i 



(14) 



(15) 



This result was first established by Pastur 14|. Zee's 
argument adds to its transparency. We also note that 
( |l5| ) is nothing but the resummed rainbow diagrams ]l5| . 
It has the suggestive form of a gap equation 



G(z) 



N 



-Tr 



z - G(z) - H L 



(16) 



familiar from mean-field (large N) approaches. In 0- 
dimensional field theories and for a Gaussian measure, 
the self-energy reduces to the two-point function. Pas- 
tur's equation is a polynomial in G of degree p+ 1. Al- 
gebraic solutions are only available for p < 3. Below we 
will show that a variety of QCD inspired random matrix 
models fulfills this restriction, with a chiral condensate 
in the massless case given by 



1 Note that the explicit form of Bd(z) is actually not needed. 
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< qq >= Im G(z — it) 



(17) 



3. The knowledge of the Blue's functions allows for a 
direct assessment of the structural character of the spec- 
tral distribution (energy bands) . This is particularly im- 
portant for systems undergoing a phase change. In gen- 
eral, a given spectrum has a support L defined by the 
sequence of branch points ±j4j, i = 1, . . . ,p with positive 
real Ai and Ax > A 2 > . . . > A p consisting of p arcs, 
[— 2j4i, — 2A2], . . . , [2A 2 , 2Ai] on which the spectrum is 
non-zero. Since the density of states at the end-point 
behaves like \J A — A, its location is determined by the 
condition 



dG 

dz 



(18) 



This condition can be translated into a condition on the 
Blue's function [[ll]. Let F(z, B(z)) be Pastur's equation 
( |l5| ) rewritten in terms of the Blue's function, i.e. 
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F(z,B(z)) = - 



1 



N ^ B{z) -a- 



-2 = 



(19) 



Then the conditions determining the end-points of the 
spectral distribution follow from the simultaneous solu- 
tions to 



dF(z, B(z)) 
dz 



= 



F(z,B(z)) = 



with 



dB{z) 
dz 



= 



(20) 



(21) 



The latter condition is just the end-point condition jl^ ) 
rewritten in terms of the Blue's function. The set of all 
end-points ±Ai in the spectral distribution comes from 
solving for B in the above set of equations ( 2p|2l| ). 

When Pastur's equation is cubic in G (Cardano class), 
the corresponding Blue's equation ( |l9|) is quadratic, and 
the end-point conditions (|2^) follow from the solutions to 
a linear and biquadratic (quartic) equation, resulting into 
four roots ±A±, ±A 2 . Therefore the spectral distribution 
consists of two arcs. The condition when A 2 vanishes 
determines a transition point, with the spectral distri- 
bution playing the role of an order parameter. Cardano 
class consists of two phases Pi and P 2 , and the spectral 
distribution is either supported by one arc (Pi -phase) or 
two disconnected arcs (P2-phase). 

When Pastur's equation is quartic in G (Ferrari class), 
the corresponding Blue's equation ([l9]) is cubic, and the 
boundary conditions reduce to the simultaneous solutions 
of a quadratic and cubic equation, thereby maximally 
six roots ±Ai, ±A 2 , ±A^. Therefore the support for the 
spectral distribution consists of either three, two or one 



arcs, corresponding to the phases P3, P 2 , or Pi. The 
critical points are determined from the condition that 
the real positive points ^3 and/or A 2 vanish. We note 
that the present analysis can be easily extended numeri- 
cally to higher degree polynomials, in the absence of the 
explicit Green's functions. This is particularly important 
for multiflavor QCD. 

4. As a direct application of the ideas discussed 
above, we now consider four different chiral random ma- 
trix models as inspired from the symmetries of the QCD 
Dirac operator and also instanton calculations. 

• Lattice with 1 flavor. 

Recently, Kalkreuter [jl(| has performed detailed lat- 
tice analysis of the distribution of eigenvalues for the 
(hermitean) Dirac operator Qi = r y 5 (p+m) using Wilson 
fcrmions in the quenched approximation. The multipli- 
cation by 75 implies that Qi = Qi. In terms of random 
matrix ensembles, Wilson fermions are represented by the 
Gaussian Unitary Ensemble (GUE) for N c > 2 and by 
the Gaussian Orthogonal Ensemble (GOE) for N c = 2. 
For one- flavor, we have 



Qi 



to 






-to 



R 



(22) 



All the entries in the deterministic matrix are N x N 
valued, while the random and hermitean matrix R is 
2N x 2N valued. Because the Wilson r-terms (needed 
to eliminate the quark doublers) break explicitly chiral 
symmetry, R is not block-off diagonal. So while the 
eigenvalues of R are not paired, the eigenvalue distribu- 
tion associated to R is symmetric about zero virtuality. 
Randomness populates evenly the positive and negative 
states. For one flavor and in the massless case, the spec- 
trum is hoped to be U(1)l x U(1)r symmetric on the 
average. The Green's function of the random matrix in 
( f22| ) for the Gaussian ensemble is given by (|J). The re- 
solvent G(z) for the deterministic and random problem 
(^||) is one of the solution to the equation Q| 



G 3 - 2zG 2 + (z 2 - m 2 + l)G - z = 



(23) 



Kalkreuter 's lattice spectra for one Wilson fermion be- 
long to the Cardano class. A typical spectral distribu- 
tion following from (|23|) is shown in Fig. [I]. It yields a 
structural change (transition from Pi to P 2 j for a critical 
mass to* = 1, with a density of states at zero virtuality 
~ A 1 / 3 Q. The localization of the heavy quark, sug- 
gests that the flavor symmetry is broken to J7y(l) with 
the appearance of a heavy quark symmetry, i.e. invari- 
ance under spin-flip. It should be clear, however, that in 
our random analysis, ( p2| ) does not account for the Ua(1) 
anomaly. The resulting spectral distribution, however, is 
qualitativaly similar to the one obtained by Kalkreuter 
on the lattice using Wilson fermions Q . 

In physical units m* = 1/E = 100-200 MeV, which is 
rather close to the strange quark mass in QCD. The ori- 
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gin of the scale S follows from the Banks-Casher relation. 
Indeed, for two flavors < qq >= — 2EA/V4 ||], where the 
chiral condensate is < qq >=<uu + dd>, with typically 
< uu >«< dd > = -((200-250) MeV) 3 . The density of 
quark zero modes is roughly n = N/V4 ~ 1 fm~ 4 . 




FIG. 1. Spectral distribution for one flavor (Wilson 
fermions) as a function of the eigenvalue A and the quark 
mass m s . 



• Finite Temperature. 

Recently, Jackson and Verbaarschot pi have sug- 
gested a simple chiral random matrix model for the mass- 
less QCD Dirac operator iJp on the cylinder, namely 



Q 2 = 



( 






R\ 


I TFT 




-A 


I 



(24) 



where irT is short for one of the two lowest Matsubara 
frequencies. Again all entries in (E3) are N x N matrix 
valued. Although the low temperature phase involves 
many Matsubara frequencies, (|24j) may be viewed as a 
schematic truncation with the right zero and high tem- 
perature content 0. 

This model belongs to the Cardano class. Indeed, 
the Green's function for the deterministic (thermal part) 
could be easily rewritten in a diagonal form by squaring 
the operator, with the result 



Gt(z) = 



{■KTf 



(25) 



Adding the corresponding Blue's function we recover the 
cubic equation ( p3| ) with the substitution m — > nT. This 
result is expected from dimensional reduction arguments 
|fL7| . At high temperature, and after suitable chiral ro- 
tations, each Matsubara mode carries a thermal mass, 
the lowest being m(T) = \/m 2 + ir 2 T 2 . The correspond- 
ing spectral distribution has been discussed numerically 
ijlql , and analytically ]lq| . The critical temperature is 
T« = 1/tt. In dimensionfull units it is T» = 30 - 70 MeV. 
Of course, this is rather unreasonable, but so is (|24|). The 
addition of more Matsubara modes affects quantitatively 
T* 0. 

The critical temperature may also be assessed analyt- 
ically by analyzing the end-points of the spectrum. In- 



deed, in the Cardano class the end-points ±A\,±A2 for 
the thermal problem are located at 



Ai = 



1 (4tt 2 T 2 - 1 + VSn 2 T 2 + If/ 2 



2 3 / 2 ttT 



V8tt 2 T 2 + 1 - I 



. l ( 47r 2 T 2 - I - V87r 2 r 2 + I) 3 / 2 

A 2 = ^/o nn 7= 5=5 I — (26) 



2 3 / 2 ttT 



^/8tt 2 T 2 + 1 + 1 



The same equations govern the evolution of the end 
points in the one-flavor case discussed above and in [Q, 
with the identification m = 7rT, as well as the end-points 
of the energy bands in the case of spin-dependent scat- 
tering off impurities in quantum Hall fluids . The real 
part of Ai vanishes at the critical point T* = 1/7T. 

• Lattice with 2 + 1 flavors. 

The issue of chiral symmetry breaking in QCD is par- 
ticularly interesting for the case of two light (up and 
down) and one heavy (strange) flavors. On the lattice 
this is usually difficult to achieve with staggered fermions. 
It is best approached using Wilson fermions. Consider 
the case where the up and down quark mass are small 
but the strange quark mass is kept fixed (2+1 flavors). 
Now, let us consider with Kalkreuter the spectral distri- 
bution for the hermitean Dirac operator Q3 = 75 {If)+ M) 
with M = diag(m„, TOd, m s ) (with m u = rrid ~ 0) and 
Wilson fermions. 




FIG. 2. Spectral distribution for 2+1 flavors (Wilson 
fermions), as a function of the eigenvalue A and the strange 
quark mass m s . 



In this case, the Green's function is given by 



G(z) = ( — Tr 



1 



1 



6N z- Q 3 



with 



Q 



3 — 



i—u,d,s 









-rrii 



R 



(27) 



(28) 



Each entry in the deterministic matrix is again N x N 
valued, while the random matrix is hermitean and 2N x 
2A^ valued. Due to the block-diagonal structure of Q3, 
the resolvent ( |2~7j ) is the sum of three resolvents, each of 
the Cardano type. Specifically, 



4 



G - -G r , 



3 G " 



(29) 



where G m are solutions to (|23| ) with the corresponding 
masses. Figure || shows the corresponding spectral func- 
tion for two light quarks m u = = 0.1, as a function 
of the heavy quark mass m s . 

For small values of m s , the spectrum is U(3)l x £^(3)_r 
symmetric, with soft breaking by the up, down and 
strange quark mass. (The assumption being that the 
r-terms do not affect considerably the chiral structure 
in the Wilson formulation). For large values of m s , the 
spectrum is U(2)l x Ur(2) with an additional invariance 
under a spin-flip of the strange quark (heavy-quark sym- 
metry). We recall that the critical quark mass for the 
Cardano class is m* = 1 or 100 — 200 MeV in physical 
units. For m s > m*, the spectral distribution resembles 
the P3 phase expected from a quartic equation (Ferrari) . 
However, the two are distinguishable by the evolution of 
their end-points. 

• Instanton Gas-Liquid 

In the instanton model of the QCD vacuum, the in- 
teraction between the instantons can cause them to clus- 
ter or dissociate into a gas, depending on the density of 
instantons per unit volume. Early simulations using in- 
stantons have revealed an admixture of free instantons 
(gas) and instanton clusters at low instanton density^. 
At a critical density n* ~ 1 fm -4 a random phase (inter- 
mediate between liquid and gas) is formed, that breaks 
spontaneously chiral symmetry p9| . Similar ideas have 
also been used to mock up the chiral transition in the 
instanton vacuum at finite temperature, using the Eu- 
clidean temporal direction as a mean of polarization p0[ . 
For these problems, the generic structure of the Dirac 
operator truncated to the space of zero modes is 



Q 



4 = 



(30) 



to d yields a Green's function for the molecules of the 
form 



1 z 
G mol (z) = a- + (1 - a ) z 2_ d 2 



(31) 



The addition law of the Blue's functions results in a quar- 
tic equation. The spectrum follows from the discontinu- 
ity of the pertinent solution to 



G = a + (1 - a)- 



G 



(z - G) 2 - d 2 



(32) 



The problem is of the Ferrari class. The roots of Q3S 
are known algebraically, so they can be used to derive 
explicitly the spectral distribution as shown in Fig. |3| 
The end-point singularities of (|32|) are given by 



Ai = ±'- 



Xj{l-d 2 )-ad 2 



(33) 



where Xi (i = 1,2,3) are the real and positive roots of 
the cubic equation 



- xT{\ + 2d 2 ) + xd 2 (d z + 3a - 1) - ad 4 = (34) 

For a = 2/3, Figure || shows a structural change in the 
spectral distribution for a hopping strength d* « or 
in physical units = 170 - 340 MeV. 



a=23 




The matrix D carries off-diagonal short (molecule) or 
medium (cluster) range order. It is N x N valued. The 
matrix elements of D in instanton models are identified 
with the hopping matrix elements between instantons 
and antiinstantons, hence chirality odd. For simplicity, 
we will restrict our discussion to the case where D is di- 
agonal since a molecule only involves one instanton and 
one anti-instanton at a time. The elements of the matrix 
D, however, need not be equal. 

Equation (^0|) is the sum of a deterministic D and a 
random R matrix. The case where D = diag(d, d, d) 
yields to the Cardano class discussed above. The case 
where some fraction a of the elements along the diagonal 
of D are zero, and all the others are non-zero and equal 



We note that this is the same as the density of quark zero 
modes, because of the Atiyah-Singer index theorem. 



FIG. 3. Spectral distribution for the instanton model 
(gas-liquid) as a function of the eigenvalue A and the hop- 
ping strength d, with fixed a = 0.66. 




FIG. 4. Behavior of the hopping strength d as a function 
of a. 
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The critical value d* depends on the concentration of 
molecules (1 — a) as shown in Fig. ||. The structural 
change occurs away from zero virtuality. 

In the vicinity of zero, (132) reduces to 



G 4 (0) + {l-d 2 )G 2 (Q)-ad 2 = (35) 
leading to the condensate (in dimensionless units) 



| < qq > | = -j=\Jl -<P + y/(l - d 2 ) 2 + Aad 2 (36) 

The sign of the condensate depends on the way we ap- 
proach zero virtuality in the complex z-plane. This re- 
sult is similar to the one discussed by Wettig, Schafer 
and Weidenmuller El]] using different arguments^. 

The critical points of (132) follow from the behavior of 
the end-points given by (|33j3^). In terms of the hopping 
strength d, the spectral distributions shown in Figs. ||, 
display novel phase structures as a function of the concen- 
tration a. The upper figure corresponds to the case where 
d is greater than the maximal critical value = For 
a = 1 we are in the totally random phase with Wigner's 
semi-circle. There are no molecules. This is the Pi phase. 
For < a < 1 we have cohabitation phases (P3 case) be- 
tween a chirally symmetric molecular phase and a spon- 
taneously broken chiral phase with a non- vanishing chiral 
condensate. For a — all the elements along the diagonal 
of D are equal to d. In this case the spectral function re- 
duces to the Cardano class, see Fig. |. The off-diagonal 
short range order through the molecules (matrix D) is 
able to destroy the long range coherence produced by 
the random distribution of R (gas). We are in the P2 
phase. The middle of Fig. || shows the spectral distribu- 
tion for the lowest critical value d, = 1. The lower of 
Fig. H shows the case with d < d* = 1. In the latter, 
the spectral distribution is always in the Pi phase. The 
off-diagonal short range order is not strong enough to 
destroy the coherence at zero virtuality. In this regime 
chiral symmetry is always spontaneously broken. 

We note that instanton models with an average instan- 
ton density of one instanton per fm 4 , display spectral dis- 
tributions that are close to the critical distribution with 
a small Jl9| ], The instanton-to- molecule model for the 
finite temperature transition in QCD [ po| , corresponds 
to the situation where a interpolates between 1 and 0, 
discontinuously in the case of a first order transition or 
smoothly in the case of a second or higher order transi- 
tion. 



3 The parameter a used in [^J corresponds to our (1 
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FIG. 5. Spectral distributions for the instanton model as 
a function of the eigenvalue A, the concentration a and the 
hopping strength d = 2 (upper), d = 1 (middle), and d — 0.5 
(lower). 



5. We have shown how the method of Blue's func- 
tions can be applied to a large class of chiral random 
matrix models. In doing so, we have unfolded two alge- 
braic equations that allow a generic classification of these 
spectra. Kalkreuter's spectral distribution for one flavor 
and Wilson fermions, as well as the massless Dirac oper- 
ator in QCD on the cylinder at high temperature (one- 
Matsubara) , have been shown to fall into the Cardano 
class (cubic equation). Wilson fermions for 2+1 flavors 
yield a spectrum that follows from a a linear superpo- 
sition of Cardano's solutions. Finally, instanton models 
with off-diagonal short range order fall into the Ferrari 
class (quartic equation). 

In terms of Blue's functions, we have discussed the 
role played by the end-point singularities and shown how 



G 



they could be used to assess the critical parameters. We 
have suggested that the spectral density of states at the 
merging end-points (zero-virtuality for instance) could be 
used as a pertinent order parameter for monitoring phase 
changes in QCD inspired spectra. This is particularly 
relevant in the massive case, where the quark condensate 
ceases to be (strictly speaking) a good order parameter. 
In doing so, we have unraveled a rich phase structure. 

The method we have followed in this paper is very pow- 
erful and allows for a number of future investigations. In 
particular, it allows for a simple assessment of the ther- 
mal aspects of the QCD phase diagram for 2+1 flavors 
using chiral random matrix models. Detailed lattice sim- 
ulations are by now available to allow for a quantita- 
tive understanding in terms of the spectral distributions. 
This issue will be discussed next. 
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